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Abstract 

Dynamical correlation functions for temperature T = are calculated for 

a Hubbard chain with infinite on-site repulsion. This chain contains three 

sites per unit cell and has a known resonating-valence-bond ground state 

for a filling of 2 particles per unit cell. A finite system of 24 sites is studied 

numerically. The recursion method is applied and results for spectral densities 

are compared to variationally calculated excited-state eigenvalues. Charge 

and spin degrees of freedom located at the backbone sites are considered. 

The spectral densities show dispersionless excitations. The energy spectrum 

displays a gap indicating an insulating ground state. Spectral functions for 

the propagation of a single hole in a resonating-valence-bond background are 

provided and indicate delocalized dispersive hole excitations. 
71.27.+a, 75.40.Gb 
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I. INTRODUCTION 

The Hubbard mo describes electronic interactions in tight binding systems. Although 
the model itself is very simple, the access to its mathematical properties turns out to be 
very difficult. The present work is based on an exactly solvable case of the Hubbard model 
in the strongly interacting limit (U = oo). For that case it has been possible to construct 
the ground state and the ground-state energy.! The class of models which are solvable by 
this method has been generalized by Strack,i Bares and Lee@ and Tasaki.iJl 

The ground state of those models exhibits the structure of a resonating- valence-bond 
(RVB) state.i One-dimensional models of this classi allowed for the calculation of equal- 
time correlation functions^ using a transfer-matrix technique. 

All studies conducted so far have been concerned with equilibrium properties of the 
ground state. Although the explicit knowledge of a Hamiltonian for the fermionic realization 
of an RVB-system enables one to study the excitation spectrum and dynamical quantities, 
no analytical statements about those quantities have been possible up to now. 

This paper presents a numerical treatment of excitations in a one-dimensional repre- 
sentative of the class of exactly solvable models introduced in [ 0. Dynamical correlation 
functions are calculated by a recursive algorithm. The resulting spectral properties are 
corroborated by a variational calculation of excited states and their energies. The method 
used here allows the calculation of single-hole spectral-functions as well. In principle these 
spectral functions are related to inverse photoemission spectra. 

II. THE MODEL HAMILTONIAN 

A linear chain with N unit cells and 3N sites is studied. The topology of this chain is 
illustrated in Fig. [l|; it may be interpreted as a chain of connected tetrahedra. The case 
of infinite on-site repulsion is assumed, which implies that only particle numbers below 
half-filling exhibit non-trivial dynamics. The ground state can be constructed for particle 
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densities of two or more per unit cell,@ thus the chain must have at least three sites (labeled 
by indices a, (3) per unit cell. 

The dimensionless Hamiltonian (the hopping matrix element defines the energy scale 
only and is set to unity) of the system under investigation is written as: 

H = p {- EE*U<w C 1 ) 

i,cr 

where c| )Q , j(J creates an electron in the unit cell i G {1, . . . , N} at site a G {1, 2, 3} with spin 
° £ {T> I}- The operators n i)0l>a . denote the corresponding occupation number operators. 
Periodic boundary conditions are assumed. As mentioned above, an infinitely strong on-site 
interaction is included by projecting all states onto the subspace without doubly occupied 
sites: 

p = IKi-^ww) ( 2 ) 

i,a 

The last term in (|l|) is a local field that lowers the energy of the connection sites a — 1. It 
ensures an almost uniform distribution of particles in the chain: The ground-state probability 
for finding an electron with spin a on a site with a = 1 (a = 2, 3) in a particular unit cell 
has been calculated as 0.3612 (0.3194) for a chain with six unit cells. 

III. THE GROUND STATE 

The ground state is derived by the same arguments as presented in [0. Consider a linear 
combination of the creation operators c\ aa : 

:= c i-1,1,ct + E C l",o- (3) 
a 

By some elementary algebra, the Hamiltonian (|T]) can be exactly converted to the following 
form 
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H = P {2N e -8N} + Y: > ( 4 ) 

i.cr 

where N e denotes the total electron number operator of the system. The Hamiltonian 
now consists of a trivial term P(2N e — 8iV) and a remainder. This remainder is positive 
semidefmite, therefore the ground-state energy is bounded from below by Eq > 2N e — 8N, 
where N e denotes the eigenvalue of the electron- number operator N e . Additionally, an upper 
bound can be derived: Any expectation value of the Hamiltonian (f|) is greater than or equal 
to the ground-state energy. The operator identity 

P*J, a P*! i(T = (5) 

allows to construct states which are eigenstates of the remainder operator ^i,a^^l,o 
with eigenvalue zero. Any state 

|*o>:=PlI*I,T*lilx> ( 6 ) 

i 

has this property. If \x) denotes the vacuum state |0), |$o) is simply a Gutzwiller-projected 
Slater determinant with two electronic orbitals in each unit cell. 

The solvability of the model is based on Eq. (|5|), a relation that holds because a local 
summation over fermionic operators vanishes due to their phase and the topology of the unit 
cell. This local property will be visible in the spectral results of the ground-state dynamics 
studied below. 

It is mentioned in passing that the norm of the states is not zero, a fact that is 
not trivial and has been proven.! For any finite system, the ground state with \x) = |0) is 
non-degenerate. This has been conjectured by Brandt and Giesekusi and proven by Tasakii 
(and Bares and Lee! for a special case), whereas a ground state with more than 2N particles 
is obviously degenerate. 

Tasakii pointed out that the ground state is generated by linear combinations of creation 
operators that contain only terms like c| jet o . — cj .cj . These terms create spin singlets 
on bonds connecting sites i, a and j, (3. Linear combinations with this structure are often 
called (short range) resonating valence-bond (RVB) states.B 



IV. SYMMETRIES 



For a numerical study of the linear chain it is essential to take advantage of symmetries. 
They allow to represent the states in terms of a significantly smaller basis set, because only 
one configuration need be stored as a representative of all states that may be generated from 
it by symmetry operations. 



The Hamiltonian is invariant under interchange of sites a = 2 and a = 3 (see Fig. 1) in 
every unit cell i. As a consequence there exists an abelian group of 2 N symmetry operations 
that commute with the Hamilton operator. All 2^ combinations of site interchanges are 
generated by every single of the 2 N following operators: 



where Mj denotes site interchange in unit cell % and I the unit operator. (The quantity e« 
can take the values ±1. Note that the product of two such operators vanishes unless they 
coincide in all signs 

It is completely equivalent to interpret the interchange of sites in terms of an interchange 
of particles, namely those particles that are located on the two participating sites. Then 
the properties of Mj acting on fermionic states become obvious immediately: Let the states 
1 00), |cr0), | Oct), \cra') denote empty, singly and doubly occupied bonds. The operator Mj 
has the following properties: 



A. Site interchange 



m({q}) = ^I1( i + *) 



(7) 



1=1 



• Mj|00) = |00) 



• Mj|cr0) = \0a), Mj|0cr) = |cr0> 



• Mi\aa') 
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From these properties, it follows immediately that the operator |(I + Mj) not only generates 
the site interchange, it additionally may be interpreted as a projector onto spin singlets for 
doubly occupied bonds. The operator |(I — Mj) projects onto triplet states. 

The symmetry of every eigenstate of the Hamiltonian may be classified according to 2 N 
different sets of signs {e^} [see Eq. (J7j)]. The ground state may be written as 

|$o) = M S |$ ). (8) 

where M s (the index s indicates that all 2 — 3-bonds are projected on the singlet-subspace) 
is an abbreviation for M({e; = +1}). This is possible, because the ground state consists 
of singlets only. However, it is important to keep in mind that excited states may not 
necessarily be written in this way. The state |$o) requires much less computer memory, 
because it is only necessary to store one representative basis state for the whole class of 
about 2 N basis states that are connected by a simple interchange of sites a = 2, 3. 

B. Translational invariance 

The second kind of symmetries that is taken explicitely into account is the translational 
invariance of the Hamiltonian. Periodic boundary conditions are assumed. Let T; be the 
operator that shifts all electrons on the chain I unit cells to the right. Again, fermion states 
are totally antisymmetric under interchange of particles, hence translation of states (wich 
may be described in terms of particle interchange) may cause a phase. For the case of I = 1 
(shift by one unit cell) the resulting sign is given by (-l)"^-™), where N e is the total 
number of electrons in the system and n the number of electrons in the last unit cell that 
are transferred to the first unit cell (i.e. moved across the "border"). The operator 

1 N 

T ? = t=E^ (9) 

VJV l=l 

acts as a generator of all translations with symmetry q = 27rm/N (m = 0, . . . , N — 1). It 
commutes with the Hamiltonian and with M s . (This is not valid for all site- interchange 
symmetries M({ej}). The operator T q obeys the relation T 9 T 9 / = S q ^T q . 
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The ground state |<3>o) may now be coded by 



$ > = T, =0 M S |$ ), 



(10) 



where |$o) contains only those basis states that may not be generated from each other by a 
translation. 

It is mentioned in passing that these symmetries lead to selection rules: The scalar 
product of states can be easily treated if one state has the site- interchange symmetry M s : 
The scalar product of the states |/) = T g M s |/) and \g) = T g /M({ej})|p) does not vanish 
if q — q' and all e; = 1 because M s and T q commute. This special case is sufficient for 
the algorithms used below. More general selection rules would have to be derived by group 
theory because the total group of symmetry operations (site-interchange and translations) 
is not abelian. 

The net storage reduction factor is of the order of N2 N (in fact, the number is slightly 
smaller, because some basis components may be eigenstates of Mj [e. g. empty bonds] or Tj) 
which reduces the multiplicity of the representative. With these considerations, system sizes 
of 24 Hubbard sites (eight unit cells) can be treated easily on a workstation. 



This section briefly recalls how dynamical correlation functions at zero temperature are 
accessible by the recursion method. A detailed discussion can be found in a very recent 
book by Viswanath and Muller.EH Correlation functions of the following type are considered: 



The operator A decribes the degree of freedom under consideration (A(£) denotes the Heisen- 
berg picture). Below, data are presented for three different A-operators, e.g. A = S? a=1 , 
which measures the z-component of the spin located at the connection sites a = 1 
with a translational symmetry q. Charge degrees of freedom are studied by e. g. A = 



V. DYNAMICAL CORRELATION FUNCTIONS 



S(t) 



($o|At(t)A(0)|$ ) 
($o|At(0)A(0)|$ ) ' 



(11) 
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exp{ig/} rii >a ^. Because the recursion method is not restricted to hermitean A-operators, 
the propagation of a single hole may be studied by the operator A = J2i expjig/} c/ ajCr . Then 
S(t) [Eq. (|TT|)1 becomes a single-particle Green function (for t > 0). 

All A-operators which are considered here commute with the symmetry generator M s 
of Eq. (|8|) because they only act on backbone orbitals (degrees of freedom regarding the 
remaining sites would have to be symmetrized such that they commute with M s as well). 
Therefore the following calculation can be performed without loss of generality in the sub- 
space corresponding to this symmetry. This can be seen immediately by formally expanding 
the state ($0!^ in terms of eigenstates with coefficients a u . Then the numerator of Eq. 
(|TT|) becomes J2v a^exp(— iE u t / h) (Q U \A\& ) where E„ denotes the eigenvalue. The matrix 
element vanishes if the symmetries of A|$ ) an d \& u ) differ. 

The method described below requires a positive semidefmite Hamiltonian. Therfore a 
shifted Hamiltonian H = H — Eq is considered where the smallest eigenvalue equals zero. For 
the hole propagator, the Hamiltonian should be positive semidefmite for particle numbers 
N e = 2N — 1. Unfortunately the ground-state energy is not known for that case, however, 
the operator H = H + AN + 2 is larger than A > 0, where A is the smallest eigenvalue of the 
non-trivial remainder in (^) for N e = 2N — 1. 

Dynamical correlation functions may be derived by an expansion of the state \^{t)) = 
A(— 1)|$ ) which solves the time dependent Schrodinger equation. This state \^(t)) may be 
expanded in terms of an orthogonal basis^ \^(t)) = J2k Dk{t)\fk)- Then the time evolution 
of the coefficients Dk(t) is governed by a set of coupled differential equations which may be 
solved by a Laplace transform. 

If this orthogonal basis \fk) is obtained by the following algorithm (which is closely 
related to the Lanczos methodQ) 



1/0) 


= A|$ ) 


(12a) 


1/1) 


= H|/o)-ao|/o>, 


(12b) 


1/2) 


= H|/ 1 )-a 1 |/ 1 >-&?|/o>, 


(12c) 
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|/* +1 >=H|/ fc )-a fc |/ fc )-6||/ fc _ 1 ), (12d) 

where 

\Jk\jk) \Jk-i\jk-i) 
the coefficient D (t) is identical to S(t) and the series of coefficients Ofc, 6| become continued 
fraction coefficients of the Laplace transform of S'(i): 

/■oo 

<*„(£):=/ dte*'S(t) 

JO 

= - • < 14 » 

4 - a 



£ - ai 



C - a 2 • • • 

If all coefficients a*,, 6| are known, the spectral density S(u) = Jf^ di exp(iotf) S(t) corre- 
sponding to the dynamical correlation function (|Tl|) can be recovered by 



S(uj) = \im 2 ft[d (u + i5)], (15) 

however, only a finite number of continued fraction coefficients can be calculated. Neverthe- 
less, the first few coefficients contain valuable information about the dynamics of the system 
which will be extracted (as will be discussed in Sec. [VTD by more sophisticated means then 
just cutting off the continued fraction and smearing out by a finite imaginary smear-out 
parameter 8. A cut-off procedure simply approximates a branch cut by a finite number of 
isolated singularities, where the number of peaks is determined by the order of the approx- 
imation. Additionally, the smear-out parameter 5 is often of the same order of magnitude 
as the hopping matrix element, otherwise reasonably looking results cannot be obtained. In 
fact, 5 must be larger than a typical distance between the isolated singularities in order to 
produce a sufficient smear-out effect. 

Finite Fermi systems have a discrete energy spectrum which only becomes continuous in 
the thermodynamic limit. Thus a recursive study of a finite Fermi system results in principle 
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in a finite continued fraction: The recursion will terminate, if the states \fk) span the Hilbert 
space completely. In practice, this termination is only seen for very small systems, because 
the number of eigenstates grows exponentially with the size of the system. 

A. Remarks on the numerical implementation 

The sequence of coefficients a^, b\ have to be computed numerically. In order to give a 
rough estimate of the numerical effort required one should consider the maximum dimension 
D of the Hilbert space: The system considered contains 3iV sites. N particles with spin up 
may be distributed on these sites in ( 3 ^) possible ways, the remaining iV spin-down particles 
may be located at the 2N unoccupied sites. Then the total dimension will be (3A)!/(iV!) 3 
which amounts to ~ 10 10 for N = 8. 

In a first step, the ground state is constructed in terms of a basis with binary coding of 
the orbital occupation. For the implementation of the recursion it is essential to generate 
the ground state with maximum precision. Here, the ground state is known, e. g. we know 
an exact algorithm which generates the ground state by successive creation of particles out 
of the vacuum, according to Eq. (g). The coefficients of the ground-state basis elements are 
coded as integers. This eliminates all rounding errors. 

The main part of the recursion Eqs. (0) is the implementation of the Hamilton operator. 
Whenever H acts on a component of a state |/), H|/) = \f), it is necessary to check 
whether the resulting basis components have already been generated before, in which case 
only amplitudes have to be added. Otherwise the newly generated basis components of \ f) 
have to be stored along with their amplitudes. This simple operation requires a large share of 
the numerical effort. In our algorithm this task is accelerated by application of binary trees. 
Then the overall numerical effort is of order D lnD', where D is the dimension of the input 
state \ f) and D' the dimension of |/'). As stated above, D and D' increase exponentially 
with system size. By taking advantage of the symmetries mentioned above, a chain length 
of eight unit cells can be treated on a workstation with 128 Mbyte of computer memory. In 
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that case D is reduced to an amount of ~ 10 6 if symmetries are taken into account. Without 
explicit use of symmetries, N = 6 is the maximum system size. 



VI. FINITE SIZE EFFECTS AND 
RECONSTRUCTION OF SPECTRAL DENSITIES 

This section provides a brief recall how the information contained in the first continued- 
fraction coefficients can be reliably retrieved. Two problems need to be addressed: First, 
because of limited memory and CPU-power only finite systems can be studied, and, second, 
only a finite number of continued-fraction coefficients can be calculated. 

The method applied here has been developed in the framework of an equivalent formula- 
tion (called Liouvillian representation) of the recursive algorithm outlined in Eqs. (0). This 
representation is based on an orthogonal expansion^ of the operator A(t) and leads to a 
continued fraction of the form 



/■OO 

CoO) := / dt e~ zt S(t) 
Jo 



1 



z H 



(16) 



a A2 

z + 



The functions S(lu) and Cq{z) are connected by a Hilbert transform: 



, \ 1 f°° , SV) 

CO * =7^ dUJ — 17 

27T1 J-oo — IZ 

S(u) =lim2 3f?[c (e-icj)] (18) 



The series of can be uniquely obtained from the afc-6|-series. 



A. Finite size effects 

Finite size effects can be eliminated by investigation of various data sets of A-coefficients 
with varying system size. The upper graph of Fig. |2| shows some data sets for chain lengths 
from four up to eight unit cells. The physical quantity under consideration is the spin 
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correlation function with A = S* =1 . The data clearly illustrate that the finite size effects 
grow with increasing index k of A*.. Therefore it is essential to restrict the analysis to 
coefficients that are in good agreement with those corresponding to shorter systems and 
ignore the size dependent higher coefficients. 



B. Termination 

The full spectral information S(u) is only accessible if all coefficients A& (of an infinite 
system) are known. The question of how much information is contained in a finite A^-series 
is by no means trivial and the goal of constructing a good approximation to the infinite 
continued fraction requires some care. 

Let Tk{z) be a terminator function that approximates the unknown rest of the continued 
fraction in "the best possible way" . 

c (z) = ^ (19) 



z + 



z + 



A 2 



z + 



z + A K T K (z) 

The crudest approximation, T K (z) = 0, is equivalent to the "cut and smear out" proce- 
dure mentioned above. More adequate terminator functions may be obtained by analyzing 
the implicit information contained in the asymptotic behavior of the Aj,-series. 

In order to construct a more adequate terminator function, the pattern of A-coefficients 



for many generic situations have been analyzed in [ 13]; the corresponding terminator func- 
tions employed here are also discussed in [ the reader is referred to these sources for 
more technical information. 

Generally, linear growth of A& with k signals a Gaussian high-frequency decay of S{ui). 
Linear growth of both A 2 k and A 2 fc+i with different slopes occurs if a spectral gap at low 
frequencies is present in S(uj). An example of this is shown in the lower panel of Fig. |^ [line 
a)]. Line b) in the same panel illustrates the characteristic change in the roles of even and 
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odd A's when an additional zero-frequency peak is present in S(u). 

The size-independent numerically-calculated continued-fraction coefficients of the system 
under investigation show the generic patterns decribed above. Thus a model spectral func- 
tion which corresponds to these patterns, S(u) = (2^/tt/ljo)Q(uj — uj 9 ) exp{ — (u — u g ) 2 /uq}, 
is expanded into a model continued fraction. The gap-parameter u g and the width of the 
shifted half-Gaussian uu are chosen such that the coefficients of the model continued fraction 
match the coefficients of the real system as close as possible. The terminator function Tk(z) 
may be recursively obtained from the Hilbert transform of ^(u;).^ 

The spectral density is reconstructed using all size-independent coefficients calculated 
numerically and - instead of just truncating - extrapolating the asymptotic pattern to infi- 
nite order by "implanting" the customized terminator function. The more size independent 
explicitly known coefficients are taken, the less important are the details of the terminator. 
The main advantage of this termination procedure is the fact that a branch cut (in this for- 
mulation located on the imaginary axis of the complex plane) of Cq(z) is not approximated 
by a finite set of isolated singularities; instead the function cq(z) displays a branch cut itself. 



The reconstruction of spectral densities is numerically performed by Eq. (0). The 
parameter e is chosen as e = 10~ 3 . Reconstruction of the spectral densities with e = 
or e = 10~ 5 yields results which coincide within the linewidth of Figs. 



VII. EXCITED STATES 

The spectral densities obtained for spin- and charge degrees of freedom indicate that 
the excitation spectrum of the chain with 2N particles exhibits a gap or an energy region, 
where the spectral densities have very little weight. In order to gain complementary infor- 
mation on the excitation spectrum of the system under consideration, excited states and 
their corresponding eigenvalues have been calculated by an iterative algorithm often called 
conjugate- gradient method.00 If the starting vector overlaps with the ground state, the 
algorithm will iteratively determine the ground state. The same procedure will yield excited 
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states and excitation energies, if it is carried out in a subspace orthogonal to the ground state 
(provided the ground state is non-degenerate). Additionally, excited states are iterativley 
calculated in different symmetry subspaces. Fig. |5] provides the energy difference between 
the ground state and the lowest excited states as a function of system size. 



VIII. RESULTS AND DISCUSSION 

The present work investigates a Hubbard model of which the ground state and its energy 
can be calculated for a special number of particles N e = 2N. This ground state has an RVB- 
structure. As far as the author knows, no other class of fermionic Hamiltonians with such a 
state beeing the exact ground state is known so far. 

However, not many of the physical properties of this ground state and the excitations 
are known. Static correlation functions [the denominator in Eq. ([Tl|)1 have been calculated 
for other one-dimensional representatives of this class! by a transfer-matrix method. In all 
cases considered so far, the static correlation functions show an exponential decay. 

The method presented above allows to calculate wave-vector dependent dynamical corre- 
lation functions making explicit use of the fact that the ground state is known. If the energy 
spectrum of the Hamiltonian has a gap, the spectral densities will show a gap as well. (The 
gap might be larger than the gap of the energy spectrum. In that case, this larger gap will 
be called dynamical gap.) The position of peaks in the reconstructed spectral densities as 
shown below reflect the dispersion relation uj pea k(q) of that particular excitation. 

In order to corroborate the conclusions about a gap in the excitation spectrum, energies 
of the lowest excited states (without symmetry restriction) and those excited states that 
couple to the ground state are calculated. 

Additionally, an attempt is made to gain information on the particle number regime 
below N e = 2N which has not been studied so far. 

Fig. |3| shows spectral densities S(u) = /f^ dt exp(iujt)S(t) [see Eq. ([Tip] for a spin 



degree of freedom A = J2i ex p(i<?0( n /,a=i,T — n «,a=i,|) which measures the ^-component of 
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the spin of the backbone electrons. All energies uo are given in units of the hopping matrix 
element (h is set to unity). The translational symmetry is characterized by q = 27rm/N 
with m = 0. . .N/2. The data indicate that there is a (dynamical) gap in the spectrum. 
There is one dominating spin excitation present with q-independent energy (dispersionless). 
Both facts support the conjecture, that the system has an insulating^ ground state and that 
the properties are dominated by local effects. Only a very small spectral feature appears for 
q — > 7r in the energy regime u ~ 3. All curves are normalized to tt. 

The operator A = J2i,a ex P(^) n i,a=i,a describes charge degrees of freedom of the back- 
bone sites. The corresponding spectral density exhibits more than just one dominating 
excitation. For q = 0, the data in Fig. |] show one very strong peak at u = 0. This simply 
reflects the q = 0-symmetry of the ground state. For other g-values, features at higher 
energies appear. Their position does not change much with q, the corresponding excitations 
seem to be rather dispersionless. 

The dynamically relevant gap is determined by the energy of the lowest excited state 
which couples via A to the ground state. Because the A-operators considered here commute 
with the singlet generator M s [see Eq. (^j)], the lowest accessible excited state must have 
the symmetry T g M s which is the symmetry of the state A(g)|$o)- The data sets 2 and 3 in 
Fig. |5| show these excitation energies as a function of the system size. They seem to saturate 
at a constant level which indicates the existence of a dynamical gap. This conclusion is 
fully compatible to the reconstructed spectral data. However, the energies of the lowest 
excited states without any restrictions on their symmetry are located within the dynamical 
gap. Data set 1 provides these energies for chain lengths up to N = 6. A final judgement, 
whether these energies saturate or not cannot be made, because larger systems cannot be 
treated without taking advantage of symmetries. 

Further, the variationally calculated energies show only very little dependency on the 
wave vector q which confirms the very small shift of the peak in Fig. [| 

The recursion method is not restricted to hermitean operators A only. Fig. ^] provides 
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results for A = J2i ex P(i<?O c ; a=i a which annihilates backbone electrons in a translationally 
invariant way. Then the quantity Sit) of Eq. ( [TTj) becomes a single-hole propagator. Its 
Fourier transform may be related to inverse photoemission spectra. 

Characteristic features of the single-hole spectral function shown in Fig. |6| are a broad 
continuum of states and a g-dependent gap. Unfortunately the origin of the u;-axis is not 
known, because the smallest eigenvalue A of the Hamiltonian H = H + 4N + 2 cannot be 
calculated. The single-hole states are delocalized (i. e. show dispersion) and may contribute 
to electronic conduction. 
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FIGURES 

FIG. 1. Two illustrations of the Hubbard chain: a) The filled circles denote the sites, the solid 
lines indicate the bonds, where hopping is possible. The dashed line frames one unit cell. Sites 
within a unit cell are labeled by an index a = 1,2,3. b) The topology of the chain may also be 
visualized by a chain of connected tetrahedra. 

FIG. 2. Continued-fraction coefficients verus k. The upper graph shows coefficients ob- 
tained for the q = spin degree of freedom A = J2i Sf a =i f° r various system sizes N = 4 — 8 with 
particle numbers of 2N. The lower graph shows two generic patterns of coefficients. Data a) are 
identical to the spin data for 8 unit cells of the upper graph. The set b) displays results for the 
q = charge degree of freedom A = J2i,a n i i a - Lines are intended to guide the eye. 

FIG. 3. Reconstructed spectral densities (Gaussian terminator, e = 1CT 3 ) for spin degrees of 
freedom A q = J2i expjigZ} Sf a=l . The first 16 continued-fraction coefficients calculated for the 
chain with N = 8 were taken into account. 

FIG. 4. Reconstructed spectral densities (Gaussian terminator [with gap for q > 0], e = 10~ 3 ) 
for charge degrees of freedom A q = J2i a exp{k//} n l a=l a . The first 16 continued-fraction coeffi- 
cients calculated for the chain with N = 8 were taken into account. 

FIG. 5. Energy differences E — Eq between excited-state eigenvalues and the ground-state 
energy for systems with N unit cells and 2N particles. Data set 1 provides results without any 
symmetry restrictions, results in set 2 (set 3) are based on the assumption that the corresponding 
excited state |^) has the symmetry Tg^Msl^) (Tq =7r M s |^')). Lines are intended to guide the 
eye. 

FIG. 6. Reconstructed spectral densities (Gaussian terminator with gap, e = 10~ 3 ) for the 
hole propagator A q = J2i exp{ig/} c x a=1 a . The first 16 continued-fraction coefficients calculated 
for the chain with N = 8 were taken into account. The origin of the w-axis is shifted by the 
unknown smallest eigenvalue A of H for N e = 2N — 1 
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